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Abstract. There has been strong observational evidence suggesting a causal con- 
nection between the binary history of neutron stars and the evolution of their magnetic 
field. In this article we discuss one of the plausible mechanisms proposed for the evo- 
lution of the surface magnetic field, that of the diamagnetic screening of the field by 
accreted material. 
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1. Introduction 



Radio pulsar observations have established that neutron stars associated with binaries have mag- 
netic fields significantly smaller than in isolated neutron stars which have field strengths clustered 
around 10 12 G (excluding magnetars). It is understood that processing in binaries alter the mag- 
netic field strength of neutron stars giving rise to the low-field binary and millisecond pulsars. 
Unfortunately, till date there is no clear consensus regarding either the nature of the internal 
configuration or the time evolution of the magnetic field in neutron stars. Depending on the gen- 
eration mechanism the field could either be supported by the crustal currents or by the Abrikosov 
fluxoids of the proton superconductor in the stellar core. 

Accordingly, two classes of models have been proposed for the evolution of the magnetic field 
in accreting neutron stars - one relating the magnetic field evolution to the spin evolution of the 
star assuming the field to be contained in the superconducting fluxoids and the other attributing 
the field evolution to direct effects of mass accretion on the crustal currents. In an accretion 
heated crust, the decay takes place principally as a result of rapid dissipation of currents due to 
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the decrease in the electrical conductivity and hence a reduction in the ohmic dissipation time 
scale. Interestingly, the mechanism of ohmic decay, unique to the crustal currents, is also used in 
models where spin-down is invoked for flux expulsion, for a subsequent dissipation of such flux 
in the crust. Therefore, the theoretical efforts in modeling the plausible mechanisms have so far 
been concentrated on the ohmic dissipation of currents in the outer crust (Konar & Bhattacharya 
2001; Bhattacharya 2002). 

Another possible mechanism is to screen the surface field by the accreting material. As the 
highly conducting accreting plasma flows horizontally from the polar caps to lower latitudes, the 
magnetic field lines are dragged along with it, by virtue of flux freezing. This dragging may lead 
to the creation of additional horizontal components at the expense of vertical ones producing an 
effective screening of the dipolar surface field. Interestingly, this mechanism would not depend 
on the location of the field in the stellar interior and would be effective irrespective of the nature 
of the interior currents. Even though this was suggested quite early on (Bisnovatyi-Kogan & 
Komberg 1974; Blandford, De Campli & Konigl 1979; Taam & van den Heuvel 1986; Romani 
1990, 1995), it is only recently that the problem is being investigated in quantitative detail. One- 
dimensional plane-parallel modeling by Cumming, Zweibel & Bildsten (2001) indicates that the 
diamagnetic screening is ineffective for field strengths above ~ 10 10 G and for accretion rates 
below ~ 1 % of the local Eddington rate. Recently, Melatos & Phinney (2001) have calculated 
the hydromagnetic structure of a neutron star accreting symmetrically at both magnetic poles as a 
function of the accreted mass. According to their calculation the magnetic dipole moment scales 
as 3 M° Ai M l a 3 where Bq, M and M a are the initial field strength, the rate of accretion and the 
total accreted mass. 



In a recent work we have proposed a 2-dimensional flow pattern of the accreted material to 
demonstrate the mechanism of diamagnetic screening (Choudhuri & Konar 2002). The accreted 
material, confined to the poles by strong magnetic stresses, accumulates in a column and sinks be- 
low the surface when the pressure of the accretion column exceeds the magnetic pressure. From 
the bottom of the accretion column, the material from both the poles move to the lower latitudes 
in an equator-ward flow, meet at the equator and submerge, pushing against the solid interior and 
displacing it very slowly in a counter-flow down-wards as well as to higher latitudes. Finally, in 
the very deep layers the material moves radially in-wards due to an overall compression of the 
star. In the top layer r m < r < r s (r s - stellar radius) of equator-ward flow, we have : 



2. Material Flow : 2-Dimensional Model 




(l) 



(2) 
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Figure 1. Flow velocity (marked a), pv, and its divergence (marked b) in a right-angular slice (0 < 9 < n/2, 
0.25 < r < 1.0) with r m = 0.75, r b = 0.5. The panels marked 1, 2 correspond to c = 0.0, 1.0. 



In the layer beneath, r/, < r < r m , where the flow turns around, we have : 



pv° r = K 2 e 
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Finally, when r < r^, the velocity is radially inward (characteristic of the radial compression in 
the deeper layers): 

pv° r =~, and pv» = 0. (5) 

It is evident from the context that the rate of accretion is related to these coefficients by the 
relation = M/4n. The coefficients K\, K 2 and K3 are related due to the fact that pv r has to 
be continuous across r — r m and r = r/,. Notice that whereas y defines the size of the polar cap 
(the angular extent being given by ~ y~ 1/2 ), B determines the magnitude of the material flow that 
sinks inward below the equator. The new material enters our region of interest only at the polar 
cap and the divergence in the top layer (r > r m ) is : 

V.(pv°) = K X e-r^-^o LZJjS. x (l + MlL^^rfi ^Scos 6)) (6) 

r \ y B sin j 

providing for a source of material only around the polar region. V.(pv) = everywhere else. 

However, with a decrease in the magnetic field the magnetic pressure in the polar region 
reduces, allowing the inflow of material through an increasingly larger region around the pole. 
Eventually, the magnetic field becomes too small to be able to channelize the flow of material 
and the accretion becomes spherically symmetric. We represent this effect through a parameter 
< Ko < 1. The effect of the widening of the polar cap is introduced through the relation 
y = (#mi n + KqA9)~ 2 where m \ n is the angular extent of the polar cap before accretion begins. In 
order to make the flow velocity more isotropic with the widening of the polar cap, an isotropic 
part is also added to the velocity field. Therefore, the expression for the velocity field, for a 
time-dependent magnetic field strength and hence a time-dependent polar cap area, is given by : 

v = (l-tfo)v (y(c)) + /r v 1 , (7) 

where v 1 is the purely isotropic part : pv 1 = (l - exp^^ 1 )) . Relating the extent of the 
polar cap with the surface field strength at a given instant of time (Konar 1997) we find that 



sin(6> min + K Q (t) A6>) 
sin 8 m \ n 



BAt) 



B,(f = 0) 



1-2/7 

(8) 



It should be noted that this relation is valid till accretion is purely spherical and Kq equals unity. 
Fig.fl] shows the profile of the velocity field and its divergence for different values of Kq- 

3. Evolution of the Magnetic Field 

The magnetic field of the neutron star evolves according to the induction equation: 

— = V x (v x B) - —V x (-V x B) , (9) 
dt 4n cr 
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Figure 2. Evolution of the mid-latitude surface field with time. Top - For a steady flow pattern with 
77 = 0.01. The curves 1 and 2 corresponds to cases without and with magnetic buoyancy (characterized by 
a radially outward velocity in the top layer taken to be v mb = 50 here). Bottom - For a time-dependent flow 
pattern with n = 0.05 and v m b = 50.0. 

where cr is the electrical conductivity of the medium. Assuming an axisymmetric poloidal field, 
allowing us to represent the magnetic field of the form, B = V x (A(r, 9)e^j, we find that A evolves 
according to the equation: 

^ + I(v.V)(sA) = ^V 2 -1Ja, (10) 

where 77 = c 2 /4ncr and s = rsinO. Evidently, it is the poloidal component of v that affects the 
evolution of A. The evolution of the magnetic field with time is obtained by integrating equation 
([Tol l subject to the boundary conditions that the field lines from the two hemispheres should match 
smoothly at the equator, requiring dA/39 — at 9 — it/ 2 and A = at 9 — such that there is no 
singularity at the pole. 

Our calculation indicates that the surface field is screened in the short time scale of the flow 
of accreting material in the top layer if magnetic buoyancy is neglected, whereas on inclusion 
of magnetic buoyancy the screening takes place in the somewhat longer time scale of the slow 




Figure 3. Field configurations corresponding to bottom panel of figure 2, at intermediate times starting 
from the initial configuration of top left-hand panel. The panels top right, bottom left and bottom right 
panels correspond to t = 0.015, 0.05, 0.1 respectively. 



interior flows, as seen if Fig. [2a]. Since, magnetic buoyancy is expected to be important in the 
liquid surface layers of accreting neutron stars, it is the second time-scale, estimated to be of the 
order of 10 5 years, which is of real importance. Remarkably, this is comparable to the time-scale 
of accretion in massive X-ray binaries indicating that the diamagnetic screening could indeed be 
one of the viable mechanisms of field reduction. 



The most interesting conclusion of our work is the fact that the screening becomes progres- 
sively less effective with the decreasing strength of the surface field as seen if Fig. [2b] and Fig. [3]. 
As the magnetic field at the surface decreases due to screening, it can no longer channelize the 
material flow and the accretion becomes more and more spherical. Magnetic field on the sur- 
face of a spherically accreting star is not subject to dragging and subsequent burial. Therefore, 
after an initial rapid decrease, the decay slows down and the field reaches an asymptotic value 
as the accretion becomes completely spherical. It should be mentioned here that similar conclu- 
sions were drawn by Konar & Bhattacharya (1997) for the case of crustal currents undergoing 
accretion-induced ohmic dissipation where a purely spherical accretion was assumed to be oper- 
ative at all times, even though the detailed crustal micro-physics, responsible for such behaviour, 
has not been incorporated in the screening model. Therefore, further investigation is needed to be 
undertaken by combining the effects of screening and ohmic dissipation taking into account the 
detailed micro-physics of the neutron star crust (Konar 2003). 



Magnetic Field of Neutron Stars 



7 



Another important aspect that requires further investigation is the question of the evolution 
of the screened magnetic field after the cessation of accretion. The effect of magnetic buoyancy, 
causing the field to re-emerge to the surface would be significant only if the flux resides in the 
topmost liquid layers. Similarly, ohmic dissipation of the field would be most effective in the 
outermost layers by virtue of smaller conductivity. It has been seen (Konar & Bhattacharya 1999a, 
199b) that for purely spherical accretion by the time accretion ceases, for almost all realistic 
binary parameters, any initial current configuration (crustal or core) would be buried deep inside 
the core of the star preventing any further evolution. Therefore, we do not expect any effective 
re-emergence of the field. However detailed quantitative estimates are still pending. 

The details of the screening model discussed here can be found in Choudhuri & Konar (2002) 
and Konar & Choudhuri (2002, 2003). 
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